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Heralded spontaneous parametric down-conversion (SPDC) has become the mainstay for single¬ 
photon state preparation in present-day photonics experiments. Because they are heralded, in prin¬ 
ciple one knows when a single photon has been prepared. However, the heralding efficiencies in 
experimentally realistic SPDC sources are typically very low. To overcome this, multiplexing tech¬ 
niques have been proposed which employ a bank of SPDC sources in parallel, and route success¬ 
fully heralded photons to the output, thereby effectively boosting the heralding efficiency. However, 
running a large bank of independent SPDC sources is costly and requires complex switching. We 
analyse a multiplexing technique based on time-bin encoding that allows the heralding efficiency of 
just a single SPDC source to be increased. The scheme is simple and experimentally viable using 
present-day technology. We analyse the operation of the scheme in terms of experimentally realistic 
considerations, such as losses, detector inefficiency, and pump-power. 


I. INTRODUCTION 

Single-photon state preparation has numerous appli¬ 
cations in the field of quantum photonics. Most notably, 
it is an essential requirement for optical quantum infor¬ 
mation processing protocols, such as linear optics quan¬ 
tum computing [1], boson-sampling [2], and optical quan¬ 
tum metrology [3] , which require high-efficiency heralded 
single-photon sources. 

Spontaneous parametric down-conversion (SPDG) has 
become the de facto standard for generating heralded 
single-photon states, owing to the relative simplicity of 
its experimental implemention. However, single source 
SPDG suffers from being non-deterministic, and generat¬ 
ing a sequence of photons in this manner is exponentially 
inefficient. To overcome this problem previous authors 
have proposed and studied various schemes for combin¬ 
ing multiple independent SPDG sources, making use of 
both spatial [4-14] and temporal [4, 15-24] multiplexing. 

Recently a fibre-loop architecture for implementing 
boson-sampling was presented [25] , where the input state 
is time-bin-encoded, and multi-photon interference is 
implemented via dynamically-controlled beamsplitters. 
Subsequently it was shown [26] that the same architec¬ 
ture can be modified to implement universal linear optics 
quantum computing. 

The viability of loop and delay-line based architec¬ 
tures has previously been experimentally demonstrated 
in the context of multiplexed, number-resolved photo¬ 
detection [27-29], quantum random walks [30, 31], and 
quantum memory [32]. These demonstrations show that 
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temporally-encoded schemes are capable of maintain¬ 
ing quantum coherence over realistic delay lengths with 
present-day technology. 

Here we show that a simplified loop-based architecture 
can implement multiplexed single-photon state prepara¬ 
tion, using time-bin encoding. The scheme requires only 
a single SPDG source with high repetition rate, a single 
fibre-loop, a dynamically controlled switch, and a single 
photo-detector with time-resolution on the order of the 
time-bin separation. The SPDG source is operated at a 
high repetition rate to minimise the length of the fibre- 
loop, and thus losses. This protocol allows a source to 
be constructed with very high heralding efficiency and 
fidelity, limited only by loss rates in the switch and fibre- 
loop, and the detector efficiency. The experimental re¬ 
quirements to build this architecture are largely available 
today, and some of the key elements have previously been 
experimentally demonstrated. 

Using a single instance of the source, loop, and detector 
immediately confers some advantages over a multiplexed 
array of such elements as we do not have to manufacture 
each element identically to ensure that the photons gener¬ 
ated are indistinguishable. Generating indistinguishable 
photons is a stringent requirement for a source suitable 
for quantum information processing. 

The concept of the scheme is to employ a single SPDG 
source, repeatedly triggered at a high repetition rate, 
yielding a correlated pulse-train in two output spatial 
modes. The pulse-train is temporally multiplexed, such 
that a successfully heralded photon is routed to one of 
the output time-bins, which then closely approximates a 
single-photon source. 

We build on a circuit layout employed by Pittman, 
Jacobs & Franson [23], extending their results by explic¬ 
itly deriving the heralding efficiency and fidelity of the 
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prepared state, taking detector efficiency, and fibre and 
switch loss rates into account. Additionally, we describe 
the optimal switching sequence to employ, so as to max¬ 
imise fidelity. We further consider the steady state be¬ 
haviour of the device, which allows single-photon prepa¬ 
ration to take place with asymptotically close to unit 
heralding efficiency, and allows immediate extraction of a 
stored photon, without having to wait for a multiplexing 
sequence to complete (i.e. an on-demand ‘push-button’ 
source). 


II. ARCHITECTURE 



FIG. 1: Architecture for multiplexed single-photon state 
preparation, as first described by Pittman, Jacobs & Fran- 
son [23]. A single SPDC source is pumped with a coherent 
state, |o), with high repetition rate, preparing a pulse-train 
with time-bin separation r. One port of the SPDC source is 
monitored with a photo-detector (either number-resolving or 
‘bucket’, i.e. on/off), which must have time-resolution bet¬ 
ter than T so as to distinguish the time-bins. The other port 
is coupled to a fibre-loop of length r, with a dynamically- 
controlled on/off switch that can be switched within time r. 
The control sequence of the switch is chosen according to the 
measurement signatures in the detected mode. Here r]s, Tff 
and r]d are the efficiencies of the switch, fibre-loop, and detec¬ 
tor respectively. At the output we have some state p, which 
describes a pulse-train in which the last time-bin ought to 
closely approximate a single-photon Fock state. 


The full architecture for our protocol is shown in Fig. 1. 
The goal is to prepare an output pulse-train in which 
the last time-bin in a sequence closely approximates the 
single-photon Fock state, |1). The pulse-train in one mode 
is measured using time-resolved photo-detection, whilst 
the other enters a loop architecture comprising a sin¬ 
gle fibre-loop, and a dynamically controlled switch. The 
dynamic switch need only toggle between completely re¬ 
flective and completely transmissive. The switching se¬ 
quence is controlled by the measurement signatures in 
the measured mode. Specifically, the switch settings at 
each time-bin are chosen so as to route the most-recent 
successful SPDC heralding event to the last output time- 
bin. Then, in the ideal case, that time-bin will contain a 
single-photon if at least one SPDC heralding event was 
successful. 

Our scheme implements a type of multiplexed state 
preparation using far more frugal resources compared to 


spatially multiplexed schemes — requiring only a sin¬ 
gle SPDC source, a single dynamically controlled switch, 
and a single time-resolved photo-detector. The architec¬ 
ture presents a significant efficiency improvement over 
stand-alone SPDC sources without multiplexing, and the 
components are largely available today. 

Below, we explicitly take into account the effects of the 
different losses in the architecture (loss being the domi¬ 
nant experimental obstacle), as well as considering both 
number-resolved and ‘bucket’ (or on/off) detectors. The 
latter is standard in present-day experiments, whilst the 
former is presently very challenging and expensive, but 
is becoming more common. The theoretical treatment we 
present could easily be extended to incorporate other er¬ 
ror models, such as dark-counts. 


III. PHOTO-DETECTION 


In the photon-number basis, a number-resolved photo¬ 
detector may be characterised in terms of a matrix of 
conditional probabilities — the probability of measuring 
m photons, given that n photons were incident on the 
detector. For a number-resolved detector with efficiency 
rjd, these conditional probabilities are given by 

Pdet(m|n) = (1) 

The pdiQt{m\n) matrix completely characterises the op¬ 
eration of the detector in the photon-number basis, and 
can easily be modified to include other effects, such as 
dark-counts, which we do not treat here. 

In the architecture we present, only two measurement 
outcomes are of interest — a single photon is detected 
(m = 1), or no photons are detected (m = 0). In this in¬ 
stance the conditional probabilities are given by 

Pdet(0|n) = (1 - rjdT, 

PdetiMn) = rid{l - ridT~^n. ( 2 ) 


For bucket detectors, which can only distinguish be¬ 
tween zero photons and more than zero photons, the con¬ 
ditional probabilities of the two measurement outcomes 


are given by 

Pdet(no click I n) 

Pdet{click|n) 


= Pdet(0|n) 

= ( 1 -%)”, 

= 1-Pdet(0|n) 

= l-(l-r?dr. (3) 


IV. HERALDED SPONTANEOUS 
PARAMETRIC DOWN-CONVERSION 

The photon source that forms the basis of our architec¬ 
ture is an SPDC source operating with a high-repetition 
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rate (i.e. being rapidly pumped). An SPDC source com¬ 
prises a crystal with a second-order nonlinearity, which 
is pumped with a strong coherent state. It prepares the 
two-mode state 

oo 

It^src) — ^ ^ 1^; ; (d) 

n=0 

and the photon-number distribution in each mode is ther¬ 
mally distributed: 

Psrc(n) = |A „|2 = 

Here n is the average photon-number, which in present- 
day experiments is intentionally kept very low so as to 
reduce the probability of multiple-pair creation events. 

The commonly measured statistics for SPDC are 
Poissonian-distributed and this distribution is frequently 
used in the analysis of multiplexed sources [4, 9, 11, 
13, 16, 18, 20, 21], often also with the probabilities 
of generating a single or multiple photons used as fig¬ 
ures of merit. The Poisson distribution arises when the 
collected light includes many spatial and/or frequency 
modes since this is the limit of the convolution of many 
thermally distributed modes. If the output of the SPDC 
modes are carefully filtered it is possible to select a single 
mode yielding a photon number distribution with ther¬ 
mal statistics of Eq. 5 [33, 34]. 

In this work we are considering a source suitable for 
quantum information tasks so the appropriate model is 
one where a single photon is prepared in a well defined 
mode so that it is as indistinguishable from other such 
photons as possible. This will allow full non-classical in¬ 
terference in subsequent photonic circuits. The fidelity of 
the photon output with an ideal single-photon state is 
the appropriate figure of merit for these tasks. If many 
modes are collected in the output state, the result will 
be a statistical mixture and this will reduce the overall 
fidelity against the ideal. As noted in the introduction, 
that using a single apparatus to generate the photons also 
naturally helps in producing indistinguishable photons. 

The key observation about Eq. 4 is that the photon- 
numbers in the two modes are perfectly correlated. Thus 
in principle with ideal detectors, if we detect some 
photon-number in the first mode, we can guarantee the 
same photon-number in the other mode. It is this corre¬ 
lation in photon-number that makes SPDC sources high- 
quality, reliable heralded sources when combined with 
high-efficiency detectors. The main problem with these 
sources is that the heralding protocol is inherently non- 
deterministic. Thus the success probability of obtaining 
m single photons from m sources for use in a larger pro¬ 
tocol drops exponentially with m. This is the motivation 
for considering multiplexed schemes, as they can signifi¬ 
cantly boost the single-photon component at the output 
of a single device. 

The photo-detection efficiency plays a major role in 
the fidelity of the prepared state. This is because detector 


inefficiency results in higher-order photon-number events 
being recorded as single-photon events, thereby preparing 
a state which is a mixture of the desired single-photon 
term, plus additional unwanted higher-order terms [35]. 
While the detrimental effect on fidelity can be reduced by 
setting n ^ 1, and hence |Ai| ^ 1, this impacts poorly 
on the efficiency. Moreover, imperfections such as dark 
counts and after-pulsing will further degrade the source. 

Next we consider the problem of combining heralded 
single-photon state preparation with our model for an 
inefficient detector. If we take a single SPDC source, pulse 
it once, and post-select on the desired heralding outcome 
in the first mode, the reduced state in the unmeasured 
mode is given by, 

oo 

Presolved ~ A/ "E Pdet(l|«)P src (n)|n)(n|, 
n=l 
oo 

^bucket = y]^>det(click|n)psrcWk)(^^|, (6) 

n=l 

in the cases of number-resolved and bucket detectors re¬ 
spectively. The normalisations are given by 

oo _ 

K = E^^det{l|«)PsrcN = 

^ (1 + n%)2 

OO _ 

y = E^det{click|n)psrcN = ■■ T- ' 

Note that these states are perfect mixtures in the photon- 
number degree of freedom, with no coherences between 
different photon-numbers. 

The conditional probabilities of preparing n photons 
given that heralding was successful, are 

Pprep('^|l) ~ I Presolved I'^) 

~ Pdet (1 |'^)Psrc ('^)/A/ 
n{n — 

~ (1 -hn)^+i ’ 

Pprep(^lclick) = (n|pbucket|^) 

= Pdet(click|n)psrc(^)/A4 
^ - [I - + Vdfi) . . 

Eor n = 1, these are equivalent to the fidelities of single¬ 
photon state preparation, conditional upon a successful 
heralding event. Note that if we were using Poissonian 
distributed output for an SPDC this would not be the 
case. 

The total probabilities of registering a successful 
heralding event are given by, 

oo 

‘^resolved — ^ ^ Pdet (I |^)Psrc(^) = A// 
n=l 
oo 

‘^bucket ~ Pdet (click|77/)psi.c (n) = A// (8) 

n=l 








4 


for number-resolved and bucket detectors respectively. 
Note that iSbucket > ‘^resolved, siucc a buckct detector ac¬ 
cepts more outcomes than a number-resolved detector. 

When number-resolved detectors of perfect efficiency 
are employed, all of the events in S correspond to per¬ 
fect single-photon state preparation. However, for bucket 
detectors or inefficient detectors, some of these events 
will correspond to having prepared the wrong photon- 
number. 


V. TEMPORAL MULTIPLEXING 


S{t) 



5 10 15 20 


t 


To increase the efficiency of a single SPDC source we 
imagine pumping it t times in succession, preparing a 
pulse-train of t time-bins but using only a single time- 
bin as output. The total probabilities of measuring at 
least one successful heralding event are. 


*^resolved (^) 


^bucket 


(t) 


1 (1 <5resolved) 


1 - 


Van 
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1 - 

1 - (1 - iSbucket)^ 


(1 


1 - 


1 

1 ^r]dn_ 


(9) 


The heralding probabilities against the number of 
time-bins t is shown in Fig. 2 for both number-resolved 
(blue) and bucket detectors (yellow), assuming perfect 
efficiency. Note that bucket detectors yield a herald¬ 
ing probability strictly greater than that of number- 
resolved detectors. As mentioned earlier this is simply be¬ 
cause bucket detectors interpret any higher-order photon- 
number term as being a successful heralding event. 
Clearly such mistaken events would degrade the qual¬ 
ity of the prepared output state so we need to consider 
the fidelity of the output as well as the probability. We 
remark that in Eq. 8 the mean photon number n and 
detector efficiency r]d only appear as a product. Indeed, 
all expressions in the paper can rewritten so that the ef¬ 
ficiency r]d only ever appears in this combination, but for 
simplicity we leave them distinct. 

When combining the SPDC source with the fibre-loop 
architecture as shown in Fig. 1, different time-bins in the 
heralded output mode may pass through the loop differ¬ 
ent numbers of times and thus suffer different amounts 
of loss. The net transmission of a time-bin subject to I 
iterations through the loop is given by, 

Ti=vJ+^r]/. (10) 


That is, the time-bin passes through the switch once, and 
then through the loop and switch an additional time for 
every iteration of the loop. 

Clearly there are many strategies that could be em¬ 
ployed when controlling the switching sequence in the 
architecture. Most trivially, as soon as an SPDC herald¬ 
ing event succeeds, we could simply keep the associated 


FIG. 2: Heralding probability (Eq. 9) against number of time- 
bins, for both number-resolved (blue circles) and bucket de¬ 
tectors (yellow squares), where 77 ^ = 1 and n = 1. 


heralded photon stored in the fibre-loop until we reach 
the last time-bin, at which point we switch it out. How¬ 
ever, this would be sub-optimal in terms of state fidelity, 
since this photon, were it to occur at an early stage, would 
be subject to many traverses of both the fibre-loop and 
the switch. 

The optimal strategy in terms of minimising loss (and 
hence maximising fidelity), is to choose our switching se¬ 
quence so as to keep the most recent successfully heralded 
photon as our final output state. 

This is achieved as follows. Whenever an SPDC herald¬ 
ing event succeeds, we couple the corresponding state 
completely into the fibre-loop, thereby coupling out what 
was previously in memory. We then toggle the switch 
so as to keep the state in memory. The state stays in 
memory unless another SPDC heralding event succeeds 
at a later stage, at which point we couple out the state 
presently in memory, and couple in the newly prepared 
state. This is repeated for the duration of the protocol. 

To calculate the fidelities of the state in the last time- 
bin with a single photon state, we must sum over all the 
ways that only 1 photon reaches the output time-bin. If 
the most recent time-bin corresponding to a successful 
heralding event passes through I loops, and a heralding 
event corresponds to n photons in this state before the 
switch, we find 

00 

^ resolved(0 “ ^ ^ Pprep(^| 1)'^^ [1 '^l] ^5 

n=l 

^ r;(l+n7?rf)^[l + n + n(l-7?rf)(l-T;)] 

[l+n-n{l-rid)il-Ti)]^ 

00 

‘^bucket(0 ~ ^ ^ Pprep('^|click)T/[1 T/] 71, 

n=l 

^ r/(l+n7/^)[l + 2n + n^77^ + nV/(l-7/^)(2-r/)] 

(1 + nTi)‘^{l+n[{l-r]d)ri + r]d])‘^ 

Here we have modelled the loss by mixing the output 
state with vacuum on a beamsplitter with a transmission 
Note that the fidelities take into account ‘accidental’ 
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single photon creation where a higher photon number is 
created due to heralding inefficiencies but then the right 
number of photons are lost to return the output state to 
a single photon state. 

For a given measurement signature, these results allow 
us to calculate the exact fidelity associated with that sig¬ 
nature. However, in reality we wish to accept all herald¬ 
ing events that are successful, and not post-select upon 
particular heralding events, which would entirely defeat 
the purpose of multiplexing. If we accept any heralding 
event that is successful, but only keep the most recent 
respective time-bin, the average fidelity of the prepared 
state is given by an average over all possible acceptable 
heralding signatures, weighted according to their preva¬ 
lence. 

Enumerating the events in decreasing order of desir¬ 
ability we label as / = 0 the case where a photon is her¬ 
alded in the last time-bin. This has probability S given 
by the appropriate expression in Eq. 8. Note the photon 
does not go through the loop in this case. The next event 
I = 1 sees a heralded photon appear in the penultimate 
time-bin and no photon in the last (or we would have 
accepted that event) and so has probability S{l—S). La¬ 
belling the time-bins in reverse chronological order from 
zero gives a probability of 5(1 — that we will accept a 
photon heralded in time-bin I which passes through the 
loop I times. Einally we append the last possible event 
where no photon is heralded in the entire train t which 
occurs with probability (1 — Sy. This forms a probabil¬ 
ity distribution for the resulting events of the switching 
strategy given by. 


p(0 


s{i-sy o<i<t 
{i-sy i = t 


( 12 ) 


Our goal is to determine when the source will yield an 
output photon and what the quality of that photon is. 
So for calculating the average fidelities we exclude the 
case where no heralding event occurred {I = t). This is 
permissible as we have a record of when the source was 
successful through the heralding signature. Excluding the 
last event means we need to normalise the distribution 
by = 1 —(1 —5)^. Thus the conditional average 

fidelities are given by. 


^resolved 


•^bucket 


t-1 

‘^resolved ^ / 7 xr-j y 

^ ^ / j '' resolved‘^resolved] 5 

Jbucket( 0 [l - 5bucket]^- 

l-(l-5bucket)^ ^ 

(13) 


In the case of perfect efficiency through the system 
{Vs = Vf = Vd = 1)5 these fidelities and heralding proba¬ 


bilities reduce to, 

^resolved 
5resolved 


= 1 , 


(1 




bucket 


1 + n’ 


^bucket 


= 1 - 


1 


1 ^ 


1 -h n 


(14) 


In that case, the fidelities are independent of t since all 
time-bins are optimal and the multiplexing sequence is 
irrelevant. With number-resolved detectors the fidelity is 
perfect, irrespective of n. Thus, the optimal strategy is 
simply to use as many time-bins as possible so as to max¬ 
imise the heralding efficiency, as pointed out in [10]. But 
for bucket detectors the fidelity is optimised in the limit 
of n ^ 0. Of course, in this limit the heralding probabil¬ 
ity vanishes. Thus, with bucket detectors there is a trade¬ 
off between fidelity and heralding probability (one very 
familiar to experimentalists developing heralded photon 
sources). 

In Eig. 3 we show the conditional average fidelity and 
heralding probabilities for different efficiencies and the 
two detectors where the source n has been optimised for 
the highest fidelity. The figure captures the effect of sev¬ 
eral parameters on the quality of the source. There is 
an overall trade-off between heralding probability and fi¬ 
delity. As the number of time-bins increases, the herald¬ 
ing probability asymptotically approaches unity, but at 
the same time the average fidelity decreases. The im¬ 
provement in the fidelity by using a heralded detector 
over a bucket detector is lost as the number of time-bins 
increases, but the main factor affecting the quality of the 
output is the efficiency of the heralding detector. 

The trade-off between heralding probability and fi¬ 
delity can be seen even more directly if we plot one 
against the other as in Eig. 4, where we use the number of 
time-bins as the parametrisation variable. Generally the 
closer to the top-right of the plot the better the source. 

In future integrated optical systems, it is plausible that 
integrated switches and delay lines might be highly ef¬ 
ficient. In this instance the fidelity of the system will 
be limited by detector efficiency. Making the assumption 
that fibre-loops and switches have perfect efficiency, the 
fidelity of the prepared states are given by. 


^resolved — 
‘^bucket ~ 


l^rjdTi 


1 2 


1 + n 

1 + 

(1 + n)^ 


These are shown in Eig. 5. Note that, 

‘^resolved ^ ‘^bucket ■> 
^resolved ^ •^bucket 5 


(15) 


(16) 


in all regimes. Eor either detector type we observe perfect 
fidelity when n ^ 0. However, the number-resolved de- 
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FIG. 3: Fidelity and heralding probabilities for various sources 
for which all specific efficiencies are the same: rf — rfs — 
Tlf — Vd- Decreasing curves plot average conditional fideli¬ 
ties J^resoived (circlcs and solid line) and *Fbucket (triangles and 
dashed line) for reducing efficiencies from top to bottom of 
?7 = 1.0 (blue), ?7 = 0.99 (green), and 77 = 0.95 (red). Increas¬ 
ing curves plot <Sresoived (t) (solid) and <Sbucket(t) (dashed) for 
same efficiencies. The shaded bands indicate the improvement 
possible by detector type, the bands themselves indicate im¬ 
provement by increasing efficiency. The optimum n for achiev¬ 
ing the highest at t = 50 was chosen for each non-unit 
efficiency detector. These were n = 0.90 and n = 0.95 for 
the resolved detector with efficiencies of 0.99 and 0.95, and 
n = 0.14 and h = 0.34 for the bucket detector with the same 
efficiencies respectively. For the unit-efficiency curves n was 
set to 0.05 for the bucket detector (which is actually max¬ 
imised for n ^ 0 ) and 1.0 for the resolved detector which 
maximises <Sresoived instead. 


tectors also yield perfect fidelity for any n if rjd = 1. This 
is not the case when using bucket detectors, since even a 
perfectly efficient bucket detector will record higher order 
photon-number terms as containing a single photon. 


VI. BLACK BOX OPERATION 

So far we have explored using the loop architecture to 
make a non-deterministic heralded SPDC source more 
deterministic. The source is pumped t times producing a 
train of t time-bins each potentially containing heralded 
photons, and switching is used to increase the probabil¬ 
ity that one of those photons will be present in the last 
time-bin in the train. We demonstrated that this protocol 
has significant benefits in terms of heralding efficiency at 
the expense of some loss of fidelity. However the result¬ 
ing source is still non-deterministic in the sense that the 
entire train of time-bins may fail to have a photon. 

The device can also be operated as a ‘black box’, 
whereby the user waits a set time and simply takes the 
output as is, regardless of any heralding signature, and we 
can determine the quality of this resulting source. Now 
the figure of merit to use is the unconditional average fi¬ 
delity T over the probabilities in Eq. 12 that also include 


FIG. 4: Gonditional average fidelity versus heralding proba¬ 
bility, parametrised by the number of time-bins (increasing 
from left to right), using sources with n set as described in 
Fig. 3. In each band the number-resolved detectors are shown 
by circles and solid line, and the bucket detectors by trian¬ 
gles and dashed line. The different bands represent different 
efficiencies from top to bottom of 77 = 1.0 (blue), 77 = 0.99 
(green), and 77 = 0.95 (red), where all specific efficiencies are 
equal (rf — rfs — rff — rfd)- More deterministic sources lie to 
the right, higher quality sources lie to the top. For compar¬ 
ison, the star in the top-right of the plot is an ideal deter¬ 
ministic single photon source and the square below it is a 
coherent state \a) with |a|^ = 1 which is sometimes used as 
a substitute. 


the no-photon outcome: 


t-i 


^resolved — ^resolved ^ ^ resolved (j ) [1 ^^resolved]"^ 5 


i=o 

t-i 


‘^bucket — ^bucket ‘^bucket (j) [1 ^bucket]'^ - 
j=0 


In the above we’ve made use of the fact that = 0 
when no photon is heralded in the entire train. These 
expressions are equivalent to the product of Eqs. 9 and 
13: JF = 5(t)J^. 

In Eig. 6 we plot as a function of n, for fixed efficien¬ 
cies (rf = rfs = rff = rfd) and a representative total num¬ 
ber of time-bins t equal to 1, 2, 4, 8 , 16, 32 or 64. Eor the 
efficiencies that we consider, it is seen that T increases 
as t increases, and that there is an optimal n that de¬ 
creases as t increases. The shift to lower average photon- 
number as the number of time-bins increases helps to 
ensure that, for imperfect efficiencies, each heralding suc¬ 
cess corresponds to a large fidelity. Eor the relatively high 
efficiencies we consider, the interplay between the maxi¬ 
mum unconditional fidelity, n, and 77 is more complex for 
bucket detectors than for number-resolved detectors, as 
heralding events that correspond to more than one pho¬ 
ton can be balanced by losing all but one photon on the 
way to the output. Nevertheless, for large n the fidelity 
reduces to .^bucket jn for all t. 

In some situations it may be the case that the rate 
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FIG. 5: Average conditional fidelity T against detector ef¬ 
ficiency Tjd and SPDC mean photon-number n, where the 
fibre-loop and switch are assumed to have perfect efficiency, 
r]f = r]s = 1. (translucent yellow surface) Number-resolved de¬ 
tectors, (solid blue surface) bucket detectors. Note that in this 
case the fidelities are independent of the number of time-bins, 
since iterations through the loop do not degrade the stored 
state. 


that photons are required is much lower than the rep¬ 
etition rate of the protocol. Alternatively we can imag¬ 
ine a bank of such sources that are left operating for a 
certain period until a large number of simultaneous pho¬ 
tons are required. Both of these cases correspond to the 
large t limit of the source. Since the heralding probability 
asymptotes to unity with increasing t, in this limit the 
source behaves as a true ‘push-button’ source, releasing 
a photon on-demand. Figure 3 suggests that in this limit 
the fidelity of the prepared state does not degrade to zero, 
but rather asymptotes to a value that is a function of the 
various efficiencies. 

Of course, one would never actually leave the device 
running in the limit of t ^ oo, since this would contra¬ 
dict the notion of ‘on demand’ state preparation, but as 
is evident from Fig. 3, more modest run times on the 
order of t 30 (in the case of r] > 0.95) are sufficient to 
bring us very close to the asymptotic state. An alterna¬ 
tive arrangement would be to have a bank of such sources 
each primed to release a single photon on the push of a 
button. In this arrangement the duration of the transient 
is not so important as the simultaneous release and in- 
distinguishability of all the photons. 

Figure 7 shows this large t behaviour for t = 
100 against n and where all efficiencies are equal 
iv = Vs = Vf = Vd)‘ Thus, Fig. 7 gives us the fidelities 
of a prepared state when the device is left running suffi¬ 
ciently long and the state in memory is coupled out on- 
demand. Observe that for low efficiencies, the bucket de¬ 
tector fidelity (yellow) out-performs the number-resolved 
case (blue). Again, this is because multi-pair events in 
which only one photon survives are able to contribute to 


the fidelity in the bucket case. 


VII. FURTHER IMPROVEMENTS 

Finally, we note some further improvements that could 
be made to the source. In Fig. 3 and Fig. 4 we opti¬ 
mised n for the conditional average fidelity at a long¬ 
time limit. This strategy can be extended to optimising 
each pulse separately — the first pulse is optimised for 
a photon heralded in the (t — 1)^^ time-bin, through to 
the last pulse optimised for the 0^^ time-bin. This could 
be accomplished by placing an additional dynamically 
controlled beamsplitter in the pump beam. 

As an example, consider = r]f = = 0.95 and 

t = 3. For these values, we calculate that /bucket reaches 
a maximum of 0.441 for the constant n = 0.668, while for 
three different pump pulses with ni = 1.31, n 2 = 0.693, 
ns = 0.466, /bucket = 0.458. In Fig. 8 we plot the maxi¬ 
mum unconditional fidelity achievable with constant op¬ 
eration, as well as this biased operation, as the number 
of time bins increases for fixed efficiencies with bucket 
detectors. This optimisation confers a larger advantage 
to bucket detectors than number-resolving detectors. 

More significantly, aside from the black box operation 
the source carries a classical record of ‘quality’ for every 
output time-bin — we know which time-bin in the train 
gave a successful heralding event. The later in the time 
train that the photon was heralded, the less loss it has 
suffered. We can then imagine having a bank of such 
sources and applying a spatial multiplexing scheme to 
pick the ‘freshest’ photons in each time-bin. 

For two parallel sources, where we pick the source that 
produces the latest heralded time-bin as the output, it 
is possible to give the joint probability distribution un¬ 
der this strategy in analytical form. The probability that 
a photon was heralded in time-bin j and not later is 
S{l—Sy. This event has to occur in at least one source, 
so there is one joint-event where both sources produce 
such photons and two lots of joint-events where only one 
source does and the other source produces an earlier time- 
bin photon, or non at all. All in all, the probability that 
the strategy produces a photon from time-bin j is 


t-i 


s{i-sy 5(i-5r + 2 ^ s{i-sy + 2 {i-sy 

k=j-\-l 

= S{1-Syy2-S). 


(18) 


After including the joint event where both sources fail to 
produce photons, the probability distribution of events 
for the two sources is 


5 ( 1 - 5 ) 2 ^( 2 - 5 ) 0 < 

( 1 - 5 ) 2 ^' 


(19) 


The resultant joint probability distribution for m 
sources running in parallel, where we take the source 













FIG. 6 : Unconditional average fidelity, T^ as a function of n for efficiencies 77 = 1.0 (top), 77 = 0.99 (middle), and 77 = 0.95 
(bottom), with number-resolved (left) and bucket (right) detectors. In each plot, the curves from bottom to top are for an 
increasing number of time-bins in the sequence t G {1, 2,4, 8,16, 32, 64}. All efficiencies are set equal, r] = r]s = r]f = r]d 


with the last heralded photon across all m sources as the 
output, can be written in a form suitable for numerical 
computation: 


Pm ('^) 


where 


T{u,ji...jrn) II PU), 
0 J Gf J 1 • ■ .Jm. } 


( 20 ) 


T{u^ Jl • • • jm) 


1 min(ji,j2,-.-im) =u 
0 otherwise 


( 21 ) 


The distributions for using up to four parallel sources are 
shown in Fig. 9. 

Using these probability distributions it is straightfor¬ 
ward to calculate the average unconditional fidelities. As 


an illustrative example, consider a loop-source that only 
has five time-bins. The improvement possible by using up 
to four such sources in parallel is depicted in Fig. 10. 

Of course these are the maximum fidelities we could 
expect, as we have not taken into account the losses in¬ 
herent in any spatial switching and routing scheme, var¬ 
ious strategies for spatial multiplexing have been previ¬ 
ously studied [5, 6, 9, 12]. Such strategies may be readily 
adapted to improving the quality of a bank of parallel 
loop sources. 

It is interesting to note that in the bucket detector 
there are parameter regimes where combining multiple 
parallel sources fares worse than just taking the output 
of a single source. Combining sources by picking the lat¬ 
est heralded time-bin reduces the effective average loss. 
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FIG. 7: Long time limiting behaviour (with t = 100) of the un¬ 
conditional average fidelities, against n and rj — rjf — rfs — Tfd^ 
where all efficiencies are equal. The translucent yellow surface 
is ^resolved, the opaquc blue surface is ^bucket- 


A7=0.99 



FIG. 8 : Maximum unconditional average fidelity, max [^], as 
a function of the total number of time bins t, for unbiased 
(blue) and biased (yellow) operation given bucket detectors, 
(top) Tj = 0.99. (bottom) 77 = 0.95. 


The parameter regimes where this strategy is a disadvan¬ 
tage are the regimes where loss is helping improve the fi¬ 
delity. The bucket heralding detector cannot distinguish 
between numbers of photons and so lets through higher 



FIG. 9: The probability distribution of last heralded photons 
in 10 time-bins for using m sources in parallel. The time-bins 
are enumerated by the number of passes through the loop. 
The probability indicated in position 11 is the probability of 
no heralding occurring. As more sources are used in paral¬ 
lel the probability of the no heralding event is reduced and 
the probability that a photon in a low-numbered time-bin is 
emitted is increased. In this example n = 0.1 and 77 = 0.95. 


order photon states, and with sufficient loss some of these 
states are reduced to single photon states on output. 

The above observation suggests that it might be ben¬ 
eficial to admit a certain fraction of two-photon states 
when using the resolved detectors. These states would 
help mitigate the effects of loss in the switching and mem¬ 
ory. The fraction of these higher-order terms accepted 
can be optimised over where in the pulse train they oc¬ 
cur and we envision quite sophisticated switching and 
parallel-source combination strategies are possible to dy¬ 
namically increase the quality of the output. Although 
outside the scope of the this document and reserved for 
future work, all of this is possible because the loop source 
produces a quality parameter for each output in the form 
of the heralding signature. 


VIII. EXPERIMENTAL VIABILITY 

Finally, let us consider the feasibility of our design with 
current and near-term technology, in terms of key exper¬ 
imental parameters including switch and fibre loss, pulse 
spreading and switching rates. We concentrate on the 
telecommunication band at 1550nm, where fibre-based 
technology is most favourable. 

In general, single-photon multiplexing is only worth¬ 
while if the individual sources produce heralded single 
photons of near perfect purity so that only one spectral 
mode is present. (This situation corresponds to the ther¬ 
mal probability distribution used earlier in the paper). In 
current 1550 nm silicon-based photon pair sources, this 
typically translates to photons with coherence time of or¬ 
der 10 — 20ps and a spatial extent of a few mm [8, 36, 37]. 

Detector efficiency and counting rate are also critical 
and there is typically a trade-off in these parameters. 
In the 1550 nm band, commercial InGaAs APD-based 
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transmission is ~ ^ 0 . 1 , at which point there is 

marginal benefit in utilising more round-trips. With the 
required loop delay of 10 ns, the fibre length is only 2.0 m 
and loop losses are negligible (below 0 . 01 % for 10 orbits) 
[42]. Similarly the dispersion length for significant pulse 
spreading is of order 20 km, and so irrelevant. For the 
foreseeable future then, the loop performance is entirely 
dominated by the switch transmission. A figure represen¬ 
tative of these parameters is shown in Fig. 11, showing 
the tradeoff between fidelity and heralding success prob¬ 
ability, parameterised by the number of time-bins. 


’O.O 0.2 0.4 0.6 0.8 1.0 

S 

FIG. 11: Fidelity versus heralding success probability trade¬ 
off, parameterised against the number of time-bins (increasing 
from left to right). Here 77 / = 1, 77 ^ = 0.8, r]s = 0.8, n = 0.5, re¬ 
flective of the parameters discussed in Sec. VIII. (blue) bucket 
detectors, (red) number-resolved detectors. 


FIG. 10: Improvement obtained by running four loop-sources 
in parallel and selecting the photons from the latest time- 
bins in each output pulse. Solid blue surfaces are for a single 
loop source, translucent yellow surfaces are for four parallel 
sources. Top graph is using resolved detectors while bottom 
graph is bucket detectors. For all configurations only five time- 
bins where used (t = 5) in each source. 


detectors operate with around rjd = 0-3 efficiency and 
detection rates of 100 MHz or better. Superconducting 
nanowire detectors offer much lower dark-counts avoiding 
an inflated heralding probability. Recent devices are also 
showing much enhanced efficiency with values of 77^7 > 0.8 

[38] , and while the maximum detection rates were previ¬ 
ously in the MHz range, they can now exceed lOMHz 

[39] . 

To avoid pulse overlap, one might then decide (fairly 
conservatively) that successive photons should be spaced 
no closer than 100 ps. This corresponds to a pump pulse 
repetition rate of 10 GHz, somewhat higher than to¬ 
day’s pair sources, and a fibre-loop length of only 2 cm, 
which would be impractically short. Instead, we consider 
a pump pulse repetition rate 1 GHz, consistent with sev¬ 
eral very recent sources [40, 41]. Assuming mean pho¬ 
ton rates of n < 0 . 1 , this is also consistent with state- 
of-the-art ceramic waveguide switches which can oper¬ 
ate at a few MHz and have already been used in spatial 
photon multiplexing experiments [ 8 ]. Such switches cur¬ 
rently exhibit transmission of order rjs ~ 0 . 8 . Significant 
additional improvements in transmission will probably be 
gradual at best. In this regime fibre losses are minimal, 
thus 77 / ~ 1. For t = 10 orbits around the loop, the net 


Another common wavelength band of interest for pair 
generation is around 800 nm where silicon avalanche pho¬ 
todiode detectors (APDs) offer much greater efficiency of 
r]d > 0.7. For most parameters discussed above, the con¬ 
clusions are largely unchanged. The fibre attenuation is 
just 3 dB/km, for which fibre losses are negligible. How¬ 
ever, the technology for high speed waveguide switches in 
this wavelength band is less mature and employing fewer 
round-trips of the loop is favourable. 

As a more generally available scenario using more stan¬ 
dard pump sources, many mode-locked lasers run at an 
80 MHz repetition rate, where successive pulses appear 
at 12.5 ns. This is a rate more compatible with the count¬ 
ing rates supported by contemporary detectors as we dis¬ 
cuss in a moment. This pulse spacing corresponds to a 
fibre length of 2.6 m. Even then, the attenuation of t = 10 
loops is of order 0.1% at 1550 nm, and only 3% at 800 nm, 
with pulse spreading still irrelevant. 

Number-resolving detectors based on transition edge 
sensors (TES) have also now been demonstrated with 
efficiencies of r]d >0.8 at 1550 nm, but the count rates 
are limited to around 100 kHz [43]. Using such detec¬ 
tors would then necessitate a matching drop in the pump 
pulse rate to 100 kHz. The corresponding fibre loop 
length of 2.1 km then gives a transmission of 77 / = 0.91 
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and the fibre and switch efficiencies begin to be compa¬ 
rable. 


IX. CONCLUSION 

We have analysed a scheme for multiplexing SPDC 
sources using temporal encoding, to allow the construc¬ 
tion of a source with higher heralding efficiency. The 
scheme requires a single SPDC source, a dynamic switch, 
a fibre-loop, and a time-resolved photo-detector. Un¬ 
like spatial multiplexing schemes, the complexity of the 
scheme presented here does not increase with the level 
of multiplexing, making it very favourable in terms of 
experimental resource requirements. The heralding effi¬ 
ciency of the system may be made asymptotically close 
to unity when enough time-bins are used. However, do¬ 


ing so comes at the expense of fidelity, which degrades as 
the number of time-bins increases. The technologies re¬ 
quired for elementary demonstrations of this protocol are 
largely available today, and could substantially improve 
the heralding efficiency of an SPDC-based source using 
present-day technology. 
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